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A dislocation moving through a quasicrystal leaves in its wake a fault denoted phason wall. 
For a two-dimensional model quasicrystal the disregistry energy of this phason wall is studied to 
determine possible Burgers vectors of the quasicrystalline structure. Unlike periodic crystals, the 
disregistry energy is an average quantity with large fluctuations on the atomic scale. Therefore the 
dislocation core structure and mobility cannot be linked to this quantity e.g. by a Peierls-Nabarro 
model. Atomistic simulations show that dislocation motion is controlled by local obstacles inherent 
to the atomic structure of the quasicrystal. 

62.20.-x,61.44.+p,62.20Fe 



In the last five years it has become possible to grow 
thermodynamically stable quasicrystals of high struc- 
tural quality and sizes j^] which have allowed one to 
measure properties like plasticity and fracture. The data 
have initiated a wealth of studies, in which the influ- 
ence of quasiperiodicity on physical properties is being 
explored. Simultaneously, potential applications of qua- 
sicrystals have been discussed, for example as friction 
reducing, oxidation resistant, non-adhesive coatings, as 
solar converters or for hydrogen storage |^ . 

The structure of quasicrystals is also responsible for 
special mechanical properties. On this subject an im- 
pressive set of experimental results is available: i) At 
room temperature quasicrystals are hard as silicon and 
extremely brittle js). ii) At about 80% of the melting 
temperature they become plastically deformable up to 
30% and soften upon straining Q. iii) Stress relaxation 
experiments reveal unusually high values for the activa- 
tion enthalpy and for the activation volume iv) For 
AlMnPd it has been proven that the deformation is gov- 
erned by dislocation motion 

In this letter we provide atomistic insight into the me- 
chanical behavior of quasicrystals. It will turn out that 
two structural features govern the motion of dislocations: 
the faults in the wake of dislocations and cluster-like ob- 
stacles in the glide planes. 

Quasiperiodic structures can be obtained as cuts of an 
irrationally oriented n-dimensional hyperplane (physical 
space) through a d-dimensional periodic crystal , where 
d > n is the number of incommensurate length scales . 
Due to the periodicity of the hyper-crystal, dislocations 
also exist in quasicrystals ^ . The d-dimensional Burgers 
vectors carry an n-dimensional component fo" in phys- 
ical space and a (d — n)-dimensional component 6^ in 
the complementing orthogonal space. A moving disloca- 
tion in a two-dimensional quasicrystal leaves in its wake a 
phason- wall, separating two areas with a phase difference 
6"*" , which gives rise to faults in the pattern and therefore 
differences in the atomic environments compared to the 
perfect quasicrystal. Atomic jumps move these faults and 
the wall may consequently vanish by diffusive processes. 



Recently, we performed numerical simulations of plas- 
tic deformation pO|JTl[ | on a two-dimensional binary 
model quasicrystal displayed in Fig. ^ We observed that 
the phason- wall remained sharp at low temperatures. 
But close to the melting temperature it was broadened 
by strain-induced diffusion and developed into a several 
bond-lengths wide glide zone, on which secondary dislo- 
cation dipoles were forming. Since this feature may be 
the central mechanism for the strain softening, we further 
study the phason-wall here. We calculate and discuss the 
"disregistry energy" , which is the potential energy nec- 
essary to shift one half of a quasicrystal rigidly over the 
other along a glide plane (line in two dimensions). Based 
on the smallest displacement that provides a minimal 
disregistry energy, we construct a generic dislocation by 
applying a kink-like displacement field. We then inves- 
tigate the dislocation core structure, the friction force 
impeding dislocation motion, and the (Peierls) stress to 
overcome it. 

We use the same binary model quasicrystal as in previ- 
ous shear simulations ||l^,|ll| . The model originates from 
the decagonal Tiibingen tiling of the plane by acute 
and flat triangles. The vertices are occupied by large 
(A) atoms. Small (B) atoms are placed in the interior 
of the acute triangles. In Fig. |l| (left) the system is de- 
picted in the atom representation, in Fig. |l| (right) in the 
bond representation, where next neighbors of different 
atomic type are connected. The model is stabilised by 
Lennard-Jones potentials. The potential depth is taken 
as I (LJ-unit) for AA and BB bonds, and as 2 for AB 
bonds in order to avoid phase separation. The LJ unit 
of length is one AB bond. 

A striking feature of the structure are concentric rings 
of ten small and ten large atoms. Due to the large num- 
ber of energetically favorable AB-bonds these rings are 
tightly bound and resemble the cluster structure of real 
systems The simple model used here therefore re- 

produces the two most important characteristics of real 
quasicrystals: namely quasiperiodicity and clusters. The 
clusters are arranged on families of parallel lines of two 
separations which form a Fibonacci sequence. The space 
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in between the large separation is traversed by a straight 
hne (dashed hne in Fig. |l|) which does not touch any clus- 
ters. This line is characterised by a low surface en ergy 
10 and serves as an "easy" path for dislocations [p^ . 
Equivalent families of lines result from rotations by 36°. 

The unrelaxed disregistry energy along such a plane is 
displayed in Fig. ^ (in LJ units). For periodic crystals 
this energy is a periodic, approximately sinusoidal curve 
with zero fault energy at multiples of the lattice vec- 
tors. For our model quasicrystal the disregistry energy 
is quasiperiodic, because the spectral decomposition for 
the spatial frequencies shows the values 1 and the inverse 
golden mean = (-s/S — l)/2. 

There are two classes of local minima in the disregistry 
energy. Those of the first class are very deep, with val- 
ues of the order of 10~^ (vertical thin lines in Fig. ||). 
The positions of these minima are projections fe" of low 
indexed hyper-crystal lattice vectors b, whose orthogonal 
component is within the acceptance domain of the trian- 
gle tiling. Therefore they are tiling vectors and represent 
the primary Burgers vectors of our model. The small but 
finite value of the minimum is the energy of the phason 
wall for phase difference b^. 

In between there are minima of a second class. Their 
corresponding values b-^ are outside the acceptance do- 
main, and the lengths fo" do not belong to the triangle 
tiling. The energies of the minima are monotonically 
(but not linearly) correlated with \b-^\. The disregistry 
energies of the second class of minima are therefore much 
higher, of the order of 1. Elastic relaxation reduces this 
value to Eph « |. 

In principle, the x-axis is densely covered with pro- 
jections fell of lattice vectors b. Therefore, a dislocation 
always can split into others of smaller Burgers vectors 
fell. However, these necessarily have larger b-^ compo- 
nents and consequently lead to a high disregistry energy. 
Such dislocations of high phasonic component will there- 
fore be connected by a high energy phason wall which 
will prevent them from splitting. 

Since a moving dislocation necessarily creates a phason 
wall with finite phase difference b-^ and finite disregistry 
energy, it is interesting to note that such a fault will be a 
preferred region for the generation of further dislocations, 
which then reduce the disregistry energy. 

In periodic crystals a Frank-Read p5[ source produces 
a sequence of dislocations with the same Burgers vec- 
tor. In a quasicrystal such a source would accumulate a 
phase difference n • n being an integer, in the glide 
plane of the emitted dislocations. As the stacking fault 
energy rapidly increases with 6^ , the source cannot oper- 
ate continuously with a constant Burgers vector. Instead 
the dislocation should split into two with smaller fe" and 
higher 6^. Only fe" components which correspond to the 
low lying minima in Fig. ^ should actually be emitted, 
while the part of the initial Burgers vector which would 
create the high energy phason-wall should be left at the 
location of the source |l^JT^. Consequently, we predict 
that a Frank-Read source in a quasicrystal will not emit 



identical Burgers vectors but Burgers vectors of varying 
length fell. It is obvious that the source will also have to 
emit Burgers vectors with larger fe^ components. This 
result is in agreement with recent experimental observa- 
tions that the probability of dislocations with larger 
fe^ components increases with plastic strain. 

Due to the homogeneity in local surroundings, in pe- 
riodic crystals the (global) disregistry energy resembles 
the local potential for a dislocatio n [|l9| and hence has 
been successfully applied by Peierls ]2Cl| and Nabarro pl| ] 
for a model of the dislocation core structure. In qua- 
sicrystals there are many different atomic neighborhoods. 
The disregistry energy is an average over the total glide 
plane. The large local fluctuations can invalidate any rea- 
soning about the dislocation core structure and mobility 
drawn in analogy to periodic crystals. This observation 
is explored by direct atomistic simulations of dislocation 
cores. 

In physical space a dislocation with phason wall was 
constructed by imposing a kink-like displacement field 
of total displacement fej = 1.9 (see Fig. ||) on a long 
slab of the binary model. The slab had a total length 
of 300 AB interatomic spacings and an aspect ratio of 
5. The dislocation was placed on a easy glide plane in 
the center of the sample. After relaxation we imposed 
homogeneous shear strains between 1% and 3%. To fix 
the deformation on the upper and lower boundary, the 
atoms were arrested in a strip four AB interatomic spac- 
ings wide. The evolution of the system then was studied 
with a constant energy molecular dynamics procedure 
from an initial temperature of 10^^ times the melting 
temperature. We calculated the displacement field with 
respect to the defect free initial configuration parallel to 
the slip plane and analyzed snapshots of the atomic sites 
taken from the simulation. 

Fig. 1^ shows the horizontal components of the displace- 
ment field parallel to the slip plane for shear deformations 
of 1%, 2% and 3% at different timesteps. At applied 
strains of 1% and 2%, there is no motion of the disloca- 
tion as a whole. Instead, it remains sessile and splits into 
two dislocations with fe^l = 1.1 and bl^ — fe^ = 0.8, which 
corresponds to the first secondary minimum in Fig. 
An increase of the deformation from 2% to 2.5% doubles 
the separation of the two dislocations. This observation 
is a first indication that not the stacking fault energy, i.e. 
the attractive force between the cores of the partials, but 
local obstacles dominate their separation. 

An increase of the applied shear deformation to a value 
of 3% leads to dislocation motion over the total length of 
the strip (third diagram in Fig. ^) . After 5000 time steps 
the dislocation is split in two dislocations similar to the 
simulation at 2% shear. Both then proceed individually 
with varying velocity and an oscillating splitting length 
until the dislocation reaches the boundary of the system. 

Our simulation therefore suggests a value of the Peierls 
stress between 2.5% and 3% of the shear modulus G. It is 
the same order of magnitude as the value calculated for 
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BCC transition metals [|19| which show brittle cleavage 
fracture at low temperature. The Peierls stress for a pe- 
riodic triangular lattice with a similar interaction model 
is about 100 times lower |^^. Thus our calculations re- 
sult in a physically reasonable value of the Peierls stress 
for a brittle material like a quasicrystal. 

Again we conclude that local obstacles on the glide 
plane dominate the dislocation motion since both dislo- 
cations moved neither at constant velocity nor with con- 
stant separation. In the following we will illuminate the 
microscopic nature of the structure intrinsic obstacles ob- 
served in our system. 

All simulations performed in this context indicate that 
the atomic environment highlighted in Fig. ^ plays the 
role of a local obstacle for dislocation motion. Here one 
large (A) atom on the upper part is bound to five small 
(B) atoms below the glide plane. This is a highly coor- 
dinated atomic cluster, which has the largest accumula- 
tion of bond energy along the "easy" plane. In Fig. ^ 
a dislocation stops in front of such an obstacle. In all 
calculations up to strains of 2% an elastic repulsion be- 
tween the dislocation core and the atomic environment 
on the glide plane with the highest coordination was ob- 
served. At higher strains, in some of our calculations, the 
leading dislocation core destroys the obstacle described 
above and transforms it into a nonplanar configuration 
which again is sessile. The dislocation is able to cut the 
obstacle only if the strain is increased beyond the Peierls 
stress. 

If the surmounting or cutting of these intrinsic obsta- 
cles is the rate limiting step in dislocation motion, the 
size of the obstacle may also explain the high activation 
volume for dislocation motion |^ and has indeed been 
proposed to do so ]2^ ]. Furthermore, it is worth men- 
tioning that an atomic cluster which is cut by the first 
dislocation is a much weaker obstacle for a second dis- 
location. Thus the cutting of the obstacles may even 
explain the strain softening which has been observed ex- 
perimentally [Q. 

Although the binary quasicrystal model used in this 
paper is very simplistic, it reflects the most important 
features of real quasicrystals, i.e. quasiperiodicity and 
clusters and therefore permits us to reach conclusions rel- 
evant and transferable to real quasicrystals. The disreg- 
istry energy of the phason wall leads to the identification 
of true Burgers vectors whose phason wall energies are 
at least one order of magnitude lower than the relaxed 
phason wall energy of these dislocations of the second 
kind. These Burgers vectors have varying length and, 
consequently a Frank-Read source should produce more 
than one dislocation type in a quasicrystal. 

Atomistic simulations were performed to calculate the 
stress which is necessary to move a dislocation. Its min- 
imal value, the Peierls stress, is of the order of magni- 
tude of other brittle materials like BCC transition met- 
als. However, our simulations reveal the dominating role 
of highly coordinated atomic environments as structure 
intrinsic obstacles for the dislocation motion. The inter- 



action between dislocations and these obstacles is cur- 
rently being studied in more detail with more realistic 
three-dimensional models. 
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FIG. 1. Binary tiling obtained from the Tubingen triangle 
tiling by decoration. Left: Atoms are displayed as disks (atom 
representation). Right: Bonds are drawn between different 
atomic species (bond representation). The ten-pronged stars 
represent tightly bound clusters and are arranged on families 
of planes (full lines). The easy planes (broken line) run in 
between the full lines of large separation. 
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FIG. 2. Disregistry energy along the easy plane of Fig. |l] 
For each minimum the orthogonal component of the Burg- 
ers vector is indicated. Erei is the energy for the relaxed pha- 
son wall corresponding to the second minimum. Lengths and 
energies are in LJ units. 
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FIG. 3. Displacements parallel to the slip plane of a dislo- 
cation at an applied shear strain of 1%, 2% and 3% and for 
different times. 




FIG. 4. Section of a configuration obtained after 50000 
time steps at a shear strain of 2%. Both partials are framed 
by Burgers circuits. The Burgers circuit around the leading 
partial is repeated in the undisturbed material below to visu- 
alize the Burgers vector. The atoms of the environment which 
acts as an obstacle for dislocation motion are drawn as disks. 
The glide plane is marked by a line. 
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